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Abstract 

A transverse multipole expansion is derived, including the longitudinal components necessarily 
present in regions of varying magnetic field profile. It can be used for exact numerical orbit 
following through the fringe field regions of magnets whose end designs introduce no extraneous 
components, i.e. fields not required to be present by Maxwell's equations. Analytic evaluations 
of the deflections are obtained in various approximations. Mainly emphasized is a "straight-line 
approximation" , in which particle orbits are treated as straight lines through the fringe fleld regions. 
This approximation leads to a readily-evaluated figure of merit, the ratio of r.m.s. end deflection 
to nominal body deflection, that can be used to determine whether or not a fringe fleld can be 
neglected. Deflections in "critical" cases (e.g. near intersection regions) are analysed in the same 
approximation . 
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I. STRATEGY AND NOTATION 



The purpose of this paper is to derive formulas for the orbit deflections caused by the 
fringe fields of non-solenoidal accelerator magnets. The main ingredient is a multipole 
expansion for fields having arbitrary longitudinal profile and including all field components 
(and only those) required to be present by Maxwell's equations. 

Because terminology describing magnets depends on context, we define some of our terms, 
if only implicitly, by using them in this section. Most magnets in accelerators are "dipoles" , 
"quadrupoles" or other "multipoles" where, in this paper, we distinguish by quotation marks 
the common names of these magnets from the dipole, quadrupole, multipole, etc., terms ap- 
pearing in mathematical expansions of their magnetic fields. The particle orbits are paraxial, 
with sma// transverse displacements, r = (x^ + with slopes {x',y') = {dx/dz,dy/dz) 

small compared to 1 because the orbits are more or less parallel to the z-axis, which is 
the magnet centerline. The dominant magnetic field components {B^, By) are therefore 
transverse to this axis, and the currents in most accelerator magnets are therefore longitu- 
dinal. But actual magnet coils must have radial leads to return the currents and, because 
of practical considerations, they also have azimuthal currents. 

The standard multipole expansion derives entirely from longitudinal magnet currents 
(this includes the bound currents in fcrromagnets). It is only for a long magnet whose 
length L is large (for example compared to a typical radial magnetic half-aperture ri/2) 
that a single multipole term provides a good approximation to the field. Yet, as concerns 
the effect of the magnet on a particle orbit, a common idealization is the short magnet or 
thin lens approximation, in which the entire defiection caused by the magnet occurs at a 
single longitudinal position. Even more extreme than our straight line approximation is to 
treat the transverse orbit coordinates {x, y) as constant through the entire magnet, body and 
ends; the deflection (say horizontal) is proportional to a field integral of the form Ax'{x, y) ~ 
jToo ^{^1 V-i ^) d'Z, where B{x, y, z) stands for any one of B^, By, dB^/dx, dB^/dy, . . . , that is, 
either of the transverse magnetic field components, or any of their derivatives with respect 
to X and/or y. Commonly then, one defines an effective magnet length L^s ~ L such that 
B{0, 0, z) dz = i3(0, 0, 0) Leff. This length is specific to the particular multipole the 
magnet is designed to produce. In spite of the facts that the magnet must be long to 
validate the multipole approximation, yet short to validate the thin element treatment, and 
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that discontinuous magnetic fields violate the Maxwell's equations, this approximation is 
curiously accurate for most accelerator magnets. Because of this good start, it promises 
to be effective to improve upon the approximation by assuming that magnets have ideal 
multipole fields within the length L^s, but also to include "end fields" applicable in regions 
of length AL_ and AL_|_ at input and output ends. In this approximation the transverse 
magnetic fields are continuous, but their derivatives are discontinuous at both ends of the 
fringe field regions. 

In a well-designed magnet, the same multipole that is dominant in the central region is 
dominant in the end regions. But the fields in the end regions are necessarily more com- 
plicated and include longitudinal components Bz{x,y,z). Since the fields in these regions 
are, in principle, constrained only by Maxwell equations, rigorous formulas for the deflec- 
tions they cause can only be evaluated by solving differential equations appropriate for the 
detailed magnet end configuration. To obtain analytic formulas we must make some as- 
sumptions, the first of which is that the formulation is not intended to apply to "intentional 
solenoids" (because of their large azimuthal currents and longitudinal field components). 
Furthermore the only longitudinal fields included are those that are required by Maxwell's 
equation to be present in regions of varying longitudinal profile. In other words, the for- 
mulas can be expected to be accurate for "well-designed" magnets, in which the dominant 
fringe field multipolarity matches the body multipolarity. This can, in principle, be assured 
by proper shaping of pole ends and proper conformation of the magnet return currents. In 
the absence of magnetic field measurements in the end regions, this is the only practical 
assumption one can make when predicting the fringe field deflections. If the fields have been 
accurately measured or calculated, to improve on formulas given in this paper, it would be 
necessary to separate out the (presumably small) extraneous components and include their 
effects perturbatively. One cannot exclude the possibility of end geometries that introduce 
multipoles for which the extraneous fringe fields are large compared to the required fringe 
fields, either intentionally of unintentionally. The present formalism would not be directly 
applicable for such fields. 

In this paper, we derive first approximations for the deflections occurring in the end field 
regions, of the form Ax'_ ~ J^^^ ^{x^y,z)dz and Ax^ ~ B{x,y,z)dz 3^. Like 

the thin lens approximation, these formulas assume the transverse orbit displacement is 
constant through the end intervals AL_ and AL_. This is a much more valid assumption 
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than assuming constant displacement through the whole magnet if, as is usually true, the 
end regions are "short"; AL± « L. Furthermore terms proportional to transverse slopes 
x' and y' can be consistently included, in the formulas for the deflections. 

A criterion for the validity of treating the end region as short can be based on the 
inequality \(3'^^y\/S.L±l f3x^y « 1., where (3x^y and are the usual beta functions and their 
derivatives with respect to the longitudinal position z. When this is true the (fractional) 
rate of change of multipole strength 1/AL± is large compared to the (fractional) rate of 
change of lattice beta functions. 

There is often a tendency to believe that multipole contributions from opposite ends of 
a magnet cancel each other. But, since this is not universally valid, in this paper no such 
assumption will be made. 



II. 3D MULTIPOLE EXPANSION 



In this section, a multipole expansion is developed that is appropriate for performing 
the calculation just described. This expansion is applicable to magnetic fields that depend 
arbitrarily on the longitudinal coordinate z but, being a power series in the transverse coor- 
dinates X and ?/, its accuracy after truncation to an order n deteriorates at large transverse 
amplitudes. The expansion is intended to describe an arbitrary "multipole" magnet along 
with its fringe field. The formalism presented here generalizes an approach described by 



Steffen and reduces to formulas he gives in the case of "dipoles" and "quadrupoles" [1|J. 

In the current-free regions to which the beams are restricted, the magnetostatic field 
B(a;, z) can be expressed as the gradient of a scalar potential $(x, z) 

B(x, 2/, z) = V$(a;, y, ^) = —x + — y + — z , (1) 

where $ satisfies 

V^$(.,,,.)^^ + ^ + ^^0. (2) 
An appropriate expansion is 

$(x,y,.) = f;f;C.,„(.)^ , (3) 

m=0 n=0 

where the coefficients Cm,n{z) depend on the longitudinal position z 3^. 
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Substituting Eq. Q into Eq. (j2I), we get a recursion relation for the coefficients; 

where in this and subsequent formulas a superscript [/] denotes / differentiations with respect 
to z\ in this case I = 2. Now, we can evaluate the gradient of the potential and get the field 
components in the three Cartesian directions 

Bx{x, y,z) = > ^ > ^ C„,,n+i{z)—, — j- 
^ — ^ ^ — ^ n! ml 

m=0 n=0 



By{x,y,z) = J2J2<^rn+hn{z)^;y^ • (5) 

m=0 71=0 
oo oo 



n\ ml 

m=0 71=0 

',,n[Z) j T 

m ml 

m=0 n=0 



The two-index coefficients Cm,n can be expressed in terms of the usual normal and skew 
multipole coefficients which, as well as being conventional, have only one index, 



[z 

x=y=0 



(6) 

x=y=0 



We next seek a representation of the field as a function of these coefficients and their deriva- 
tives. The relation can be applied recursively to obtain 

where the upper limit of the series k is equal to the integer part of m/2. This shows that 
the coefficients Cm,n can be expressed as a series of even derivatives of Co,n+i or Ci,„. Using 
Eq. (jni) we can distinguish two cases for m, namely m = 2k (even) or m = 2/c + 1 (odd), 
and we have 

k 



1=0 ^ ^ 

C2k^l,n = J2(-^)i% 
;— n \ / 



n+2k-2l 

=0 

The requirement Co,o = corresponds to the restriction to non-solenoidal magnets. 
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Substituting this representation into Eqs. (0) and rearranging the m-summation yields 



oo oo TO 



Bx{x,y,z) 
By{x,y,z) 



n=0 m=0 

oo oo 

EE 



n=Q m=0 1=0 

oo oo 

n=0 m=0 



n! (2m)! 



/ / n\ (2m 

E 

./=o 

m+l 

E 



J2«] 



+ a 



[21] 

n+2m-2l 



^ / "n+2m-2l 



1=0 



m + l\ [21] y 
I r'^+2-+i-2'2m+l 



(9) 



n—n m— n ;— n \ / V / \ 



+ a, 



[2i+l] 

n+2m-l-2i 



n=0 m=0 1=0 

again hmiting the ranges so the lowest coefficients are bo = Ci^ and = Co,i. 

In an idealized model of a magnet, only one (or in the case of combined function magnets, 
two) of the multipole coefficients will be non-vanishing in the body of the magnet (length 
Lcs) and in this region only the / = terms in the expansions survive. The important 
terms are: (m = 0, / = 0) corresponding to the leading "design" multipole; (m = 0,/ = 1), 
the "next-to-leading" term associated with longitudinal variation of the design multipole; 
and (m = 1, / = 0) coming from the next higher body multipole. Examples in this paper 
are mainly concerned with the relative importance of the ffist two of these terms in the 
deflections caused by the actual magnet, including body and ends. The same formulas 
could, however, be used to evaluate the relative importance of the second and third terms — 
to answer the question "Which are more important, fringe fields or body field imperfection?" 

To obtain results concerning the symmetries of the skew and normal multipole coefficients 
it is more useful to express these formulas in terms of cylindrical coordinates. This is done 
in Appendix A. 

In the fringe regions of the magnet, the fields can be arranged so that they match the 
central fields at the ends of the body region and fall linearly to zero in the fringe regions. 
For example, let us keep just one more term as a "next approximation", arrange its leading 
(/ = 0) part to match a given body field at ^ = 0, and let it vary linearly with z; 



Bx{x,y,z] 
By{x,y,z] 
Bz{x,y,z] 



^n— 1 y2m+l 



oo oo 

'^^^ (n- l)!(2m + l)! 

n=l m=0 ^ ' ^ ' 

°o °° ™n „,2m 

^ ^ ^ n\ (2m 



-1^ 



r,[0] , dl] 

"^n+2m ~ "n+2m ^ 



n+2m 



+ 6^'^ z 



(10) 



n=0 m=0 

oo oo 

EE 

71=0 TO=0 



r^n y2m+l 

n\ (2m + 1)! 



.[1] 

"n+2m 



where the n index has been shifted by 1 in the expansion for convenience in the next 
step. Next, we arrange for Bx{x, y, AL) = by setting 

iM 

,[1] _ 'Jn+2m /-.-.N 
"n+2m — ■ \^^) 

It can be seen that this condition also assures By{x, y, AL) = 0. This is a consequence of 
the requirement that V x B = 0. Setting B(x, y,z) — Q for z > AL, we have assured that 
the transverse field components are continuous. Due to the artificial assumption of linear 
fall-off of the field in the fringe region, the longitudinal component Bz is discontinuous in 
this approximation. 

At this point, the "multipole" magnet has been idealized by a model whose parameters, 
apart from its multipolarity index, are its multipole strength b^n+2mi lengths Les 

and AL-I-. This representation is appropriate for representing the magnet within a parti- 
cle tracking computer program. The lengths AL-t could be determined by best-fitting to 
measured fringe fields. But, to reduce the number of parameters in the remainder of this 
paper, and with some reduction in accuracy, a slightly different approach will be taken; the 
impulses delivered by the fringe fields will be evaluated in a way that is independent of the 
fringe field lengths: all the integrals involved will be computed by using the "hard-edge" 
approximation, i.e. taking the limit for which AL± 0. In this limit the straight line 
approximation becomes exact. 

For the sake of consistency another point must also be made. Since the dominant mul- 
tipole in the magnet body is also dominant in the fringe field, there can be an appreciable 
contribution to the dominant field integral (due to the magnet as a whole) that comes from 
the fields in the fringe regions. It is a matter of taste whether this contribution is to be 
treated as part of the main field or part of the fringe field. In this paper, from here on, to 
simplify the formulas somewhat, the term "fringe field" will refer to components other than 
the dominant component, but restricted to those components necessarily associated with the 
dominant multipole. In other words, the contributions from the dominant multipole compo- 
nent in the fringe regions will be counted as part of the ideal magnet field integral. Treating 
the magnet in this way increases its effective length probably making it more nearly equal 
to the the physical magnet length; i.e. L Les, and this will be assumed in all subsequent 
formulas. 
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III. DEFLECTIONS AT MAGNET ENDS 



For a given magnet with a perfect 2(ri+l)-pole geometry written in cylindrical coordinates 
(see Appendix A), the scalar potential satisfies the following symmetry condition: 

^r,e,z) = ^r,^-e,z) , (12) 
n + 1 

which leads to a relation between the harmonic multipole number allowed by symmetry n' 
and the multipole order (n + 1): 

n' = (2j + l)(n + l)-l . (13) 

Thus, for a normal "dipole" {n = 0) the multipole coefficients allowed by the magnet sym- 
metry are of the form b2j, for a normal "quadrupole" (n = 1) &4j+i, for a normal "sextupole" 
[n = 2) &6j+2; etc. Consider now a "multipole magnet", with normal symmetry, for example. 
Following the symmetry condition (fT^ . we can rewrite the field components ()A7jl . keeping 
terms of the expansion to leading order: 

BA., y, z) =Zm { (^±M!Mf) _ (- + -^)""[(- + ^)x-^.(. + l),]6L-k^) ^ ^ 4) 

nir . A f (^ + ^yr&n(^) (x + Zy)"+M(^+l)^-^(^ + 3)y]&l''(^) , p. , 

ay\x,y,z)-Ke<^ n! 4(n + 2)' +c^l^ + 4j 

E.(x, z) =Jm I ^''^y^l)!"^ + + 

(14) 

where the functions 0(j) represent polynomial terms in the transverse variables x, y of order 
greater or equal to j. These expressions apply for n > 0. The special case of the "dipole" 
will be treated separately. Here the terms proportional to and approximate the fields 
present due to the longitudinal field profile variation and do not include fields that could be 
present due to non-ideal magnet design. 

For a particle traversing the magnet along the straight line having transverse coordinates 
{x,y), the impulse {i.e. change of transverse momentum) imparted by the nominal field 
component is 

7^e{(x + i?/)"} 



^pI = - e / By{x,y,z)dz ^ - e6„Lefr- 

ApJ = e / B^{x,y,z)dz ^ ebnL^s j 

body ^• 
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where L^g = J bn{z)dz/bn is the effective length of the magnet, and fe„ is the nominal field 

body 

coefficient in the body of the multipole magnet. The quantities in Eq. fll5j) . the intentional 
and dominant ("zero order") deflections caused by the magnet, are only approximate, since 
they account neither for orbit curvature within the body of the magnet nor for end field 
deflections. Expressions like this will be used only as "normalizing denominators" in ratios 
having (the presumably much smaller) magnet end deflections as numerators. For magnets 
other than bending magnets, for which the average deflection is zero, it will be necessary to 
use r.m.s. values for both the normalizing denominator and the numerator. 

The impulse due to the fringe field at one end of a magnet is defined in this paper as 
the effect of field deviation from nominal, from well inside (where the nominal multipole 
coefficient is assumed to be independent of z) to well outside the magnet (where all field 
components are assumed to vanish). These will be the limits for the integrals used in order to 
calculate the fringe deflection. To obtain explicit formulas the upper limit of these integrals 
will be taken to be infinity. Exploiting the assumed constancy of x and y along the orbit, 
these integrals will all be evaluated using integration by parts. 

Suppressing the entire pure multipole contribution, as explained above, we have 
f^^B{x,y, z)dz ^ 0. For x = y = this is an equality by definition, and for finite 
displacements it is approximately true if, as we are assuming, the transverse particle dis- 
placements remain approximately constant. This is consistent with our straight line orbit 
approximation. 

The individual components of the impulse can themselves be separated into terms due 
to longitudinal fields (labeled ||) and due to transverse fields (labeled ±); 

Api^^ = Api^iW) + Api^yil.) , (16) 

where 

Ap^(||)= e / y'B,{x,y,z)dz 

fringe ^-j^r,-^ 

^P^(II) = - e / x'B^{x,y,z)dz 

fringe 

are the momentum increments of the particle caused by the longitudinal component of the 
magnetic field and 

^Px(^) = - e / By{x,y,z)dz 

fringe ^^g^ 

^P^(-L) = e / B^{x,y,z)dz 

fringe 
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are the momentum increments of the particle caused by the transverse components of the 
magnetic field. Using the leading order expressions of the magnetic field, we obtain the 
relations 



eh 

Ap/(||)^ ^-^Im {{x + zyr^'jy' 



eh 



(19) 



and 



Ap^(±) ^ — ^7^e {(a; + iyY \{n + l)xx' + (n + ^)yy' + i{n - l)xy' - i{n + l)yx']} 

4(n + 1)! 



Ap^(±) ^ — — ^^——Im {{x + iy)"" [(n + 3)xx' + (n + l)yy' + i{n + l)xy' — i{n — l)yx']} 

(20) 

The total impulses caused by the fringe field are therefore 

- 4(^_['^), ^e{(^ + ^VT [{n + l)(x - iy){x' + ly') + 2iy\x + ly)]} 

■ (21) 



Ap^ ^ — — ^^^-—Xm{{x + iyY[{n + l){x — iy){x' + iy') — 2x\x + iy)]} 
^ 4(n + 1)! 

Even though they occur at a fixed point in the lattice, because these impulses depend 
on slopes x' and y' and are truncated Taylor series, they are not symplectic. To use them 
in long term, damping-free tracking, symplecticity would have to be restored by including 
deviations in transverse coordinates 



IV. APPLICATION EXAMPLES 

The formulas just derived are appropriate to calculate the end field deflection of any single 
particle. But to assess the importance of these deflections it is appropriate to calculate their 
impact on the beam as a whole, for example by calculating an r.m.s. deflection, such as 
{Ap-^j_)rms = \j ((Apl)^) + {{Apy)"^). Here the operator (.) denotes an averaging over angle 
variables. Note that here, and from here on, the subscript ± specifies the transverse impulse, 
and does not refer to a magnetic field component. Formulas for r.m.s. values like these are 
derived in Appendix B. This section contains examples of the use of those formulas, starting 
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with the cases of flat and round beams, then speciahzing the results further for "dipole" and 
"quadrupole" magnets. The derived formulas are flnally applied for evaluating the impact 
of magnets end flelds in the case of the Large Hadron Collider (LHC) and the Spallation 
Neutron Source (SNS) accumulator ring. The calculations are based on Eq. ()B13|) . 



A. Flat Beam 



For a flat beam, one of the transverse degrees of freedom {e.g. the vertical y, y') vanishes. 
Thus, the total transverse r.m.s. momentum increment from the magnet body is 




(Ap^)_ = yWpm - ' Wei , (22) 



where /3" represents the average of the (3"' in the body of the magnet and e± is the transverse 
emittance. The total transverse r.m.s. momentum increment from one of the fringes of the 
magnet is 



where f3 and a represent the beta and alpha functions, at the fringe location. The ratio of 
these quantities is 

{Ap{\^s e± /(2n+l)/5"[l + (2n + 3)a2 



(24) 



(ApD^ms SLeffY (n + l)(n + 2)/3'^ 

Assuming that the beta functions are not varying rapidly, if the magnets are in non-critical 
locations (which is to say most magnets), the square root dependence can be neglected, so an 
order-of-magnitude estimate (dropping an n-dependent numerical factor not very different 
from 1) is given by 



eff 



(25) 



The case in which fringe fleld deflections are likely to be most important is when a is 
anomalously large, for example in the vicinity of beam waists such as at the location of 
intersection points in colliding beam lattices. In this case, (again dropping a numerical 
factor) the ratio of deflections is roughly 



{Ap{ 



rms 



. ^ - • (26) 



rms -'^eff 
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The same result is obtained by setting j3x » Py in Eqs. flB13|) . 

Often the relative deflection is so small as to make neglect of the fringe fleld deflection 
entirely persuasive. The simplicity of the formula is due to the fact that the fringe contri- 
bution is expressed as a fraction of the dominant contribution. Note that, as stated before, 
this formula applies to each end separately, and does not depend on any cancellation of the 
contributions from two ends. In fact, nonlinear analysis shows that in magnets fringe-fleld 
contributions can tend to add up instead of cancelling Q]. 



B. Round beam 



For a round beam, the two transverse emittances are equal ex = ey = e±. For simplicity, 
we assume that typical values of horizontal and vertical lattice functions are approximately 
equal; Px ^ Py = (3 and ax ^ ay = a. Also assume that ~ i.e. the beta functions 
do not vary significantly in the body of the magnet. Taking into account the previous 
hypotheses, the total transverse r.m.s. momentum increment for the body becomes: 

1/2 



n/2 



5^2(1/2, -n, -n; 1, 1/2 - 1 



(2r2- 1)!! 



(27) 



where the function in the square root represents the generalized Hyper-geometric function 
(see 0| for details). Applying the same simplifications, the r.m.s. momentum kick given by 
the fringe field is: 



(Api). 



2"+3(n + 1)! 



./=o 



2{n - I) 
n — / 



1/2 



(28) 



where we considered I3x ^ (3y = (3 and the same for the a functions. Notice now that the 
sum of the coefficients gn,i = Qn^ifl + 9ni,i + 5'n,i,2 depends only on a^. The series involving 
them can be also written as a sum of a few generalized Hyper-geometric functions. The 
ratio of the r.m.s. momentum transverse kicks is: 



(Ap 



'_|_)rms 



CJa" 



(29) 



where the coefficient C„ is: 



[n + l] 



n!ELo(t-?)e/)^n, 



5^2(1/2, -n, -n; 1, 1/2 - n; l){2n - 1)!! 



1/2 



(30) 
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Let us consider two before: one where a is small and one where a is large, as near 

the interaction points of large colliders. For the first case {a small), we may neglect the 
terms having a as a factor in the coefficient gn,i and in the second case, we can pull out a 
from the square root and neglect terms in the coefficient Qn^i having now the a function in 
the denominator. In this way, the coefficients C„ of Eq. (|3(J|) will depend only on the order 
n. We plot in Figs. ^ the behavior of these coefficients as a function of the multipole order 
n, for large and small a. The dominant factor in C„ seems to be l/{n + l), which is reflected 
in the slow asymptotic decay depicted at the plots. For all practical cases (multipole orders 
up to 20), Cn lies between 1/2 and 1/10. Assuming now that the average /? in the body of 
the magnet is not so different from (3 in the fringe, one gets for small a functions: 



as in Eq. (j2Sl), and for a large: 



as in Eq. 



rms 



(Api) rms 



, (32) 
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C. Dipole magnet 

Consider a "straight" dipole magnet; the configuration of poles and coils is symmetric 
about the x = and y = planes, and the coils are excited with alternating signs and equal 
strength. By symmetry 5^ is odd in both x and y, By is even in both x and y, and B^ is 
even in x and odd in y. Using the general field expansion of Eq. Q, we get: 



°° ™ ( 2n+l 2m+l /™ , , , 

" ^ ^ (2n + l)!(2m+ 1)! V/ > 



m,n=0 1=0 
oo m 



2n+2m+2-2l 



- ^l^^l^ ' {2n)\{2m)\ ^ / J^2n+2™-2^ • (33) 

R - \^ I Ij X y /m\ [2;+!] 

2^ 2^ (2n)!(2m + 1)! V ^ 

m,n=0 1=0 ^ > ^ < J \ / 

Taking the field expansion up to leading order, we get: 



B, = b2xy + 0{4) 
5, =2/ 6^+0(3) 



By = bo-h^^V + lb,{x'-y') + 0{4) , (34) 
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where 62 represents a sextupole field component allowed by the symmetry of the "dipole" 
magnet (for an ideally designed magnet 62 = 0) and 0(3) and 0(4) contain all the allowed 
terms of higher orders. 

A point has to be made about the application of the integrals evaluating the rms momen- 
tum kicks for bending magnets: because of the curved central orbit, these integrals are not 
exact, as previously mentioned. Nevertheless, in most practical cases, the field uniformity in 
the interior of a "dipole" magnet is very high, and thus, on heuristic grounds, this approach 
can be expected to provide fairly good estimates even in this case. 

The change of transverse momentum imparted by the dipole field is (see Eq. fll5|l ) 



Ap = — e y bodz ^ —eboLcs , (35) 

body 

where as before Lcs = /body ^o'^'^/^o effective length of the "dipole" magnet, and 60 is 

the main dipole field in the body of the "dipole" magnet. Using Eq. p8|) the deflections in 
one fringe are 

Apl 2eboyy' , Ap^ ^ -eboyx' , (36) 
and the total r.m.s. fringe kick is 



{ApiUs = ebo^J A{y^y'') + {y^x") . (37) 
Using Eqs. (|B10jl and (IBlljl . we have 

ivY) = , {y^x'^) = {y^) {x'^) = + f^^y'^'y , (38) 

and the r.m.s. transverse momentum kick becomes 



Thus, the by-now-standard ratio is 



(Arf)rms ^ _l_ /(1+MM + (1 + c^DPye.ey 

{Ap\\ms^ 8 4/3, ■ ^ ^ 

Except for numerical factors near one this formula yields the same "ball-park" estimates as 
given by Eq. and Eq. for the small a and large a cases. 
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D. "Quadrupole" magnet 



The configuration of poles and coils in a "quadrupole" magnet is symmetric about the 
four planes x — 0; y — 0; x — y; x — —y and if the coils are excited with alternating signs 
and equal strength, the magnetic field will satisfy the following symmetry conditions: 
is even in x and odd in y; By is odd in x and even in y; B^ is odd in both x and y; and 
Bz{x, y, z) = Bz{y, x, z). As before, we may express the field components as: 



^-^Y. Yl (2n)!(2m + l)! \l 1^' 

m,n=0 1=0 ^ ^ ' J \ 



m 

2n+2m+l-2l 



^2/ - 2^ 2^ (2n+ l)!(2m)! V I i '^2n+2m+i-2/ ■ (41j 

m,n=0 1=0 ^ / V / \ 



m,n=0 1=0 

The field expansion can be written as 



2n+2m+l-2l 



h-^{3x' + y')b? 



(42) 



^ y bi-—{3x' + y')br +0(5) 

By = X bi - —iSy' + x')b'^^ +0(5) ' 
B, = xyb?+0{4) 

where bi{z) is the transverse field gradient at the quadrupole axis, and 0(4), 0(5) contain 
all the higher order terms. For a particle traversing the magnet with a horizontal deviation 
X and vertical deviation y from the center, the momentum increments produced by the 
nominal field gradients are 

Apl = -ebixL^s , ApJ = ebiyL^s , (43) 

where L^s — /body ^i^"^/^! effective length of the quadrupole magnet. The momentum 

increments of the particle contributed from the longitudinal component of the magnetic field 
are 

APx(ll) ~ exyy'bi , Ap^(||) ^ -exyx'bi , (44) 
and the momentum increment produced by the transverse component of the fringe fields are 



Api{±)^—^[2xyy'+ix' + y')x'] , Apf{±)^-^[2xx'y + ix' + y')y'] . (45) 
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Combining the contributions, the total momentum increments due to fringe field are 

Apl [2xyy' - {x^ + y^)x'] 

eh • ^^^^ 

Again, by averaging the sum of squares of the transverse momenta contribution, we obtain 
the total rms transverse momentum kick imparted by the fringe field: 

e6~ f 3 

{/^pilms ~ ^ I (1 + 5a^)/3,e^ + ^ [(1 + - 8a.a,/3,/3, + 2(1 + Sa^)/?^] ,2^^ 

3^ 1 ■ 

+ (1 + bal)l3yel + — [(1 + al)(3l - Sa^ayp^py + 2(1 + 3al)(3l] e,ej| 

(47) 

Note that the expected rotation symmetry of the quadrupole is exhibited both in this formula 
and in the body deflection formula. The standard ratio is 

(ApjU. 1 [ (1 + 5al)PlPyel + 3/?. [(1 + a^Pl - Sa^yMy + 2(1 + ?>al)f3l] ejey 
^ (1 + hal)Ml4 + Wy [(1 + - 9,a.ayMy + 2(1 + ?>al)l3l] e.e^ \ 

2My{fxea^ + Wy^y) J 

(48) 

Again dropping factors near 1, this leads to the same ball-park estimates of Eq. and 
Eq. dSl . 



E. Magnets of LHC and SNS 

The LHC and the SNS accumulator ring are good examples for testing the validity of 
the derived fringe fleld figure of merit formulas. Indeed, the purpose of these two proton 
machines and thereby their magnet design differs in great extent: the LHC, a high-energy 
hadron collider, is filled with long super-conducting magnets of very small aperture (around 
1 cm). In contrast, the SNS ring, a low-energy high intensity accumulator, contains short 
normal conducting magnets with wide aperture (tens of cm). In addition, the lattice design, 
optics functions and physical parameters of the two machines are substantially different, e.g. 
the emittance of the SNS beam is several orders of magnitude bigger, than the one of the 
LHC. In Table m we summarize the parameters of the main magnets in the two accelerators 
entering in the figure of merit formulas ()4U|) and (PS|) . 
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In Fig. 121 we plot in logarithmic scale the fringe-field figure of merit estimates for the 
LHC and the SNS accumulator ring magnets. The black bars represent evaluation with the 
exact formulas derived for dipoles and quadrupoles (see Eqs. and (011)) and the grey 
bars represent the evaluation with the formula for round beams (j3H) . In both cases, the 
total effect for each magnet is computed by summing up the fringe-field figures of merit from 
both ends due to all the magnets of the same type. The fringe field importance in the case 
of the SNS is striking, especially for quadrupole magnets, whereas in the case of the LHC 
can be completely neglected. Note that similar results can be derived by careful dynamical 
analysis and computation of tune-shifts due to fringe fields or dynamic aperture analysis 
for both the LHC [3] and the SNS j^. It is important to stress that even the approximate 
formula for round beams (j3ip is slightly pessimisitic and within a factor of 2 of the exact 
figure of merit. 



V. CONCLUSION 



We have derived formulas for the momentum kicks imparted by the fringe fields of general 
straight (non-solenoidal) multipole magnets. These formulas are based on an expansion 
having arbitrary dependence on the longitudinal coordinate. This expansion can be used for 
direct integration of the equations of motion for particle tracking or other analytical non- 
linear dynamics estimates. It also permits the fringe part and the body part of individual 
magnets to be identified and separated. A figure of merit, the ratio of r.m.s. end deflection 
to r.m.s. body deflection is introduced and evaluated. Its proportionality to the transverse 
emittance results in an easily-evaluated measure of the importance of fringe flelds both in 
cases in which the variation of optical functions is not too rapid and in the opposite case 
of rapid variation. These results are in agreement with previous crude estimations which 
employed simple physics arguments based on Maxwell laws jojl. Finally, the formalism has 
been applied to the most common cases of multipole magnets, namely normal "dipoles" and 
"quadrupoles" 13]. Since the straight line approximation has been used throughout, these 
formulas are only precise for magnetic flelds that are well-approximated by step functions 
(the "hard-edge" approximation). Thus, the formulas contain no parameters associated with 
the fringe shape (for example, see Ulll^l). Also, as stated previously, only those fringe flelds 



matching, and therefore required by, the nominal body multipolarity are accounted for. 
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Numerical evaluation of the end/body figure of merit shows that fringe fields can be 
neglected in the magnets populating the arcs of large colliders like the LHC. In these rings, 
the magnets are long enough and the emittances are so small (of the order of 10^^ m rad) 
that the effect of fringe fields is a tiny perturbation as compared to the dominant multipole 
errors in the body of the magnets. The effect may be important, however, in small rings, 
as the SNS accumulator ring or the muon collider ring 3|; where the emittance is large 
(typically 10~^ m rad) and the magnets much shorter. Careful consideration should be also 
taken in the case of the magnets located in the interaction regions of the collider where 
the beta variation is quite big. 

It is perhaps appropriate to call attention to possible "overly optimistic" use of the 
scaling law. Often quadrupoles are grouped in doublets or triplets in which the desired 
focal properties rely on the intentional, highly-tuned, near cancellation of deflections caused 
by more than one element. In such cases, the fringe deflections are, of course, amplified, 
when evaluated relative to the gross multiplet deflection. This effect is most obvious at focal 
points. 

Since the early analytical studies of Lee- Whiting [l^ and Forest , significant progress 
has been achieved for the construction of accurate maps which represent the motion of 
particles through the magnet fringe field, using either direct numerical evaluation with exact 



integration of the magnetic field 



16 



17| or parameter fit of an adequate function 



m 



IS 



{e.g. the Enge function j2jj). These maps are essential for the study of non-linearities 
introduced by fringe-fields through Hamiltonian perturbation theory techniques. On the 
other hand, the scaling law we have emphasized can provide a rough estimate of the impact 
of these fringe fields in a ring. If the fringe fields are found to be important, a thorough 
numerical modelling and analysis of their effect has to be undertaken, including computation 



m 



of the amplitude dependent tune-shift, resonance excitation and dynamic aperture 

innrinr, 

121 . 122L 12 3L 1241 . 125L I26| , as non-linear dynamics can be very sensitive to the details of different 
lattices and magnet designs. Furthermore, great care is required to preserve symplecticity 
and use these maps in particle tracking. 
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APPENDIX A: 3D MULTIPOLE EXPANSION, CYLINDRICAL COORDI- 
NATES 



The magnetic field representation in Cartesian coordinates {x, y, z) is not optimal for 
studying symmetries imposed by the cylindrical geometry of a perfect multipole mag- 
net. For this, it is preferable to rel y o n expansions in cylindrical coordinates (r, 9, z) = 



(a/x^ + y^, arctan (y/x), z) 



oie to rel y o n expansions 



Both expansions are equivalent and 



the use of the former or the latter depends mostly on taste and the specific problem to be 
treated. 

First, consider the magnetic scalar potential written in the following form 



oo oo 



I. n=0 m=0 ) 

where now the ^-dependent coefficients Qn+i m{.z) are generally complex. The above expan- 
.on follows erectly from the fact that the Laplacau commutes wth a/90 1^. Th,s allows 

the consideration of solutions where the dependence in 9 is an harmonic 2{n + l)-pole. This 



expansion is compatible with the general so^ 



Q,lai,Q 



ution of the Laplace equation in cylindrical 



coordinates, involving Bessel functions 

Using Eq. (jAljl and the Laplace equation, one gets that Qn+ifi = 0. Moreover, Qu+iaj 
should vanish for n > (all terms except the dipole). Finally, we have a recursion relation [3|, 
Ql similar to Eq. 0: 

Qn+i,m+2{z) = 7 — , .y^'") for m^n-l , (A2) 

[n + 1)^ — [m + ly 

where again the superscript in brackets denotes derivatives with respect to z. Following 
these relations, one can show that all coefficients with m < n + 1 vanish. Thus, the first 
non-zero coefficient is Qn+i,n+i (for m = n + 1). By extending the recursion relation ()A2|) so 
as to express any coefficient as a function of Qn+i,n+i, we get: 



[2k] 



The summation indexes can be rearranged so as to express the magnetic scalar potential in 
cylindrical coordinates 0, 0| : 




^,-i(^),n+l+2fc. ^ (A4) 
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and the three-dimensional field components are: 



B (r d z)-ne j Ve'("+^)^ V + g^^k+i] 

B,[r, ^.z)-lte<^ye ^ ^^^^^ + l + k)\k\ ^"^^ 



{z) r 



n+l+2fe 



The coefficients Qn+i = Qn+i,n+i can be related with the usual multipole coefficients, through 
Eqs. (ini). First, we write the scalar magnetic potential in Cartesian coordinates: 



{-lf{n + 1)! 

2 

.n=0 fc=0 

The magnetic field components are computed by the gradient of the potential ()A6|) : 



\. TL J — 



[{n + l + 2k)x-i{n + l)y]GZ\{z) 
[x' + y^f~\x + iyr''^x . (A7) 



[-(n + l)x + z(n + 1 + 2A;)|/] G^^^^{z 



Using Eqs. (0), we get: 



n/2 



fc=l 

"/2 / -.Nfc/ 



an{z) = in + 1)! 7^e{^?„+l(z)} + ^ 22^(71 + 1 + fc)!fc! ^e{^n+i-2fe(^)} 

(A8) 



where the upper limit of both series is the integer part of n/2. Thus, in the absence of 
longitudinal dependence of the field, the normal and skew multipole coefficients are just 
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scalar multiples of the imaginary and real part of Qn+i{z). On the other hand, the situation 
is more complicated in the case of 3D fields. By inverting the series f)A8|l . we have: 

n/2 



k=0 

n/2 



(A9) 



k=0 



-2k\^) 



^^=0 .... ^^^^^ 



where the coefficients 7?.^\ and TZ^^k '^^^ t)e computed order by order by the j + 1 relations 

T^nor ^ 1 V- {-inn + l-2k)in + l-2k)\ ^ 

(n + 1) ' ^ 22^(n + l + A;)!A;! '^n-^kj-k 

1 ^{-lf{n + l-4:k){n + l-2k)\ 

k=0 

and j runs from 1 to the integer part of n/2. Using the last relations, the scalar potential 
and the magnetic field can be expressed as a function of the usual multipole coefficients. By 
expanding the complex polynomials in the expression of the magnetic field components, one 
recovers the expansions of the magnetic fields Q in Cartesian coordinates. 



(n + 1) '4^ 22fc(n + l + A;)!fc! '^n-2k,j-k 



APPENDIX B: EVALUATION OF R.M.S. END DEFLECTIONS 



In order to evaluate the r.m.s. defiection caused by a magnet end, we start from the 
expressions by sphtting the product inside the brackets: 



4(n+ 1)! 



TZe {{x + iy)^} [(n + l)xx' + (n — l)yy'] 
+ Im {{x + iy)"'} [— (n + 3)xy' + (n + l)x'y)] 



(Bl) 



4(n+ 1)! 



TZe {{x + iy)""} [(n + l)xy' — (n + 3)x'y] 
+ Xm {{x + iyy} [{n — l)xx' + (n + l)yy') 



The total r.m.s. transverse momentum kick imparted by the fringe field is (Apj 



{{Aply) + {{Aply), where the operator (.) denotes the average over the angle variables. 
An equivalent expression stands for the defiection due to the body part of the field. The 
(.) operator is linear, we can first compute the sum of squares of the momentum kicks and 



21 



then proceed to their averaging. Thus, we have: 



(Api) 



rms 



4(n+ 1)! 



(A {ne{{x + iynf + h (Im {{x + iynf 

^1/2 

+2/37^e {{x + iyr} Im{{x + iy)"}) ' 



ebnLes 



ni 



[((7^e {{x + iynr + {Im {{x + iyT}f)] 



where /i, /2 and /s are: 

A = (n + 1)^x2 (x'' + + y2 [(^ + 3)2^/2 + _ i)2^/2j _ ^ 

/2 = x2 [(n - 1) + (n + 3)2y'=^] + (n + lfy'^{x''^ + 7/'^) - 8(n + l)xx'yy' 

/s = 4 [-(n + l)(x2 + y2)a;Y + a;^(a;'2^^'2-jj 
We have the following relations for the real and imaginary part of {x + iy)": 



[n/2] 

ne{{x + ^yr}= E(-l)\2/) 



and thus: 



[(n-l)/2] 



{ne{{x + iyr}f 



[{x' + y^ + neiix + iyr}] 
1 " 

E 



1=0 



(Xm {{x + zy)^)' = [{x' + y^)" - Tee {(x + ^i/)^"}] 



n 



1=0 



^2n-2ly2l 



ne{{x + iyY}Im{{x + iyY} = -Tm {(x + iy)^"} 



^2n-2i-1^2i+l 
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where the upper hmit of the last sum is taken to be / = n for uniformity in the equations, 
instead of the last non-zero term for which I = n — 1. Finally, it is straightforward to show 
that 



2 /™2 , „.2 



(7^e {{x + iyT}) + iTm {{x + iyY})' = (x 



l=Q ^ ^ 



(B6) 



After expanding the products in Eq. ()B2|) and collecting the terms of equal power in the 
transverse variables, we have that the transverse kicks can be written in the following form: 



where the Qh's are 



eb„ 



A{n + 1)\ 



1=0 



1/2 



.1=0 



1/2 



= {u;,{nJ)+u;2in,l)){x'--''+'x"){y'^) 



(B7) 



= (^3(n,0+^4(n,/))(x2"-V)(y2'+2) 



= Mn,l)+cu,{nJ)){x'-^''+'){yV) 
= {u;,{n,l)+u,{n,l)){x'''-'^){y'^+Y) 



(B8) 



= uj{n,l){x'^-''-'x'){y''+y) 



with the coefficients cu^'s: 



{u;r{n,l)+u;s{n,l)){x'-~'^^'x'){y''^'y') 



uji{n,l) 
uj5{n,l) 



K + 1) 



n 



(n^ + 4n + 5) 
-2 {n + 2){-iy 



2n 
21 



LJ2{n,l) 
ujQ{n,l) 



-8{n + l){n - l){-iy f2n 



2/ + 1 



21 



2n(-l)' 



2n 
21 



2(5n + 2/n + 2)(-l)' f2n 

21 + 1 V2/ 
-2/(2n + f2n 

21 + 1 V 2/ 



uJs{n,l) = -8(n + 1) 



(B9) 
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In order to proceed to the averaging of the transverse variables, we write them in the standard 
form 

{x, y} = V (^x,y(3x,v Ca:,y, {x', y'} = J^{S^^y + a,j,^yC^^y) , (BIO) 

where Cx^y are the transverse emittance associated with the corresponding phase space di- 
mension, (5x,yi OLx,y are the usual beta and alpha functions and C^, Sq stand for cos 0q, sin 0^, 
respectively. Using the above relations and averaging over the angle variables 0^ one can 
show that: 



m ) 22'" 



2m\\l + {2m+l)a'l\j3^-^e 
m J 



m+l 



22™+i(m + 1) 



(Bll) 



Then, the flkS become: 



1^1 = {u;i{n, I) + LJ2{n, /)) 



m + l / 22'"+2 



2{n - l)\ /2A (2n - 2Z + 1)[1 + (2n -21 + 3)al](3^-'(3, 
n-l )\l) 22^^+2(71- Z + + 2) 



2]f:{n-lfql n-l+2 I 
y X 



^ , , ^ n^ /^2(n - l)\ f2l\ {21 + !)[! + (2n -21 + l)c^^Pr'-'P^'^^^'^^ 



n 



22'»+2(n-/ + l)(/ + l) 



^ , , , (2{n - l)\ (2l\ {2n - 2/ + 1)[1 + {21 + l)a^]/:^r'+^/3;-^6»-^+^6j; ^ 

Q3=Mn,0+-5(n,0)( JiJ 22'^+2(,-/ + l)(/ + l) 



Q4 ^ {uji{n,l) + ujQ{n,l)) 



n — I J \ I 



2{n - l)\ (2l\ {21 + 1)[1 + {21 + 3)a2]/5^-'/3^e^-'e'+2 



22"+2(/ + !)(/ + 2) 



J J 22'^-^2(/ + + 2) 



t ( i n^ ^2(n - l)\ (2l\ {2n -21 + 1){21 + l)a.a,/3r^/3^6r'+^6i; ^ 
Mn,l)+us{n,l))^^_^ 2^n^^{n-l + l){l + l) 



(B12) 
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After collecting terms of equal emittances, the r.m.s. transverse momentum kicks can be 
expressed as: 



{Ap{) 



eb„ 



2n+3(n + 1)! 



.1=0 ^ / \ / 



m J2—m 

y 



1/2 



ebnL 



eff 



_|_ jrms 



2"n! 



1=0 



1/2 



(B13) 

where the bars on the /5's denote their average values over the body of the magnet. The 
coefficients gn,i,m, given by 

[{n' + 1){21 + 1) (;) - 2/(2n + 1)(-1)'(5)] [1 + {21 + 3)a^] 



gn,l,oiS^x,y^ 



x,y> 



9n,l,l{0ix,y, Px,y) 



(/ + !)(/ + 2) 

8(n + l)(n - /)(2/ + 3)(-l)'(',';)a^ay/3^ 
+ !)(/ + 2) 

[(n^ + 4n + 5) (2/ + 1) (;) + 2(5n + 2/n + 2)(-l)' ('^;)] [1 + (2n - 2/ + 

/3,(n-Z + l)(/ + l) 

[(^2 ^ 4^ ^ 5)(n^) _ ^ 2)(-l)^fi;)] (2n - 2/ + + (2/ + l)al](3x 
^ /3j,(n-/ + l)(/ + l) 

8(n + 1) [(2/ + l)^;) + in- a] (2r2 - 2/ + 

(n-/ + l)(/ + l) 



9n,l.2iS^x,yi Px,y) 



[{n^ + 1) (7) + 2n(-l)'(g)] (2n - 2/ + + {2n -21 + 3)ag] 
(n-/ + l)(n-/ + 2) 



(B14) 



depend on the twiss functions a^^y, (3x,y and on the multipole order n. One may note 
that r.m.s. transverse momentum kick of the fringe is represented by the square root of a 
polynomial of order n + 2 in the transverse emittances and ty as compared to the square 
root of a polynomial of order n representing the body contribution (see also P]). Thus, their 
ratio should be proportional to the transverse emittance. This scaling law is indeed exact 
for the case of the "dipole" and "quadrupole" . For higher order "multipoles" , it is exact for 
flat and round beams (Sec. IV). 
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Eq. (jnni))) wlien tlie a function is small (left) and when the a function is large (right). 
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FIG. 2: Fringe-field figure of merit estimates for the LHC and the SNS accumulator ring magnets. 
The two different bars in each case represent evaluation with the exact formulas derived for the 
fringe-field figure of merit of Eqs. (|4U|) and (|48|) (black bars) and the approximate formula for round 
beams (|lTT|l (grey bars). 



TABLE L Parameters associated with the LHC and SNS magnets, whose fringe-field figure of 
merit is evaluated in Fig. |2 When two numbers occur, they are associated to the minimum and 
maximum value. 



Magnet 


Number 


-^^eff [m] 


/3x,y 


[m] 




[m] 


\^x,y 


1 [m] 


€x,y [m rad] 


LHC Quadrupole Triplets 


16 


5.5 - 6.37 


1055 - 


4463 


1157 - 


4401 


1.1 - 


203.9 


5.03 10-1° 


LHC Arc Quadrupoles 


368 


3.1 


32 - 


178 


32 - 


176 


0.5 - 


- 2.4 


7.82 10"^ 


LHC Dipoles 


1104 


14.3 


28 - 


176 


40 - 


143 


0.5 - 


- 2.6 


7.82 10"^ 


SNS Dipoles 


32 


1.5 


4 - 


8 


6 




1.1 - 


- 1.9 


4.8 10"^ 


SNS Quadrupoles 


52 


0.5 - 0.7 


2 - 


28 


2 - 


26 


- 


- 8 


4.8 10"^ 
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